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ABSTRACT 

We perform a joint likelihood analysis of the power spectra of the 2dF Galaxy Red- 
shift Survey (2dFGRS) and the cosmic microwave background (CMB) anisotropics 
under the assumptions that the initial fluctuations were adiabatic, Gaussian and well 
described by power laws with scalar and tensor indices of n s and n t . On its own, the 
2dFGRS sets tight limits on the parameter combination f2 m ft[^J but relatively weak 
limits on the fraction of the cosmic matter density in baryons f2f,/f2 m . The CMB 
anisotropy data alone set poor constraints on the cosmological constant and Hubble 
constant because of a 'geometrical degeneracy' among parameters. Furthermore, if 
tensor modes are allowed, the CMB data allow a wide range of values for the physical 
densities in baryons and cold dark matter (uji, — fl^h 2 and u> c = fl c h 2 ). Combin- 
ing the CMB and 2dFGRS data sets helps to break both the geometrical and tensor 
mode degeneracies. The values of the parameters derived here are consistent with the 
predictions of the simplest models of inflation, with the baryon density derived from 
primordial nucleosynthesis and with direct measurements of the Hubble parameter. 
In particular, we find strong evidence for a positive cosmological constant with a ±2<r 
range of 0.65 < fl\ < 0.85, completely independently of constraints on J7a derived 
from Type la supernovae. 

Key words: Galaxy clustering, large-scale structure, cosmic microwave background- 
cosmology: miscellaneous. 



* Here h is Hubble's constant Ho in units of lOOkms 1 Mpc 1 . 
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1 INTRODUCTION 

Until recently, cosmology was a subject starved of data, with 
poor or non-existent constraints on fundamental quantities 
such as the curvature of the Universe, the power spectrum 
of density irregularities, and the cosmic densities in baryons, 
cold dark matter and vacuum energy. The situation has 
changed dramatically over the last few years. Following the 
discovery of the CMB anisotropics (Smoot et al. 1992) it 
was realized that many of the fundamental parameters of 
our Universe could be determined via accurate, high resolu- 
tion measurments of the CMB (e.g. Bond et al. 1994, Jung- 
man et al. 1996). This has now become a reality through a 
number of exquisite ground based and balloon experiments 
(see Halverson et al. 2001; Lee et al. 2001; Netterfield et 
al. 2001). Constraints on cosmological parameters derived 
from these experiments are described in several recent pa- 
pers (de Bernadis et al. 2001; Pryke et al. 2001; Stompor et 
al. 2001; Wang, Tegmark & Zaldarriaga 2001). 

Significant advances have also been made in survey- 
ing large scale structure in the Universe. The development 
of wide-field correctors and multi-fibre spectroscopy means 
that it is now possible to measure redshifts of hundreds of 
thousands of galaxies. Two such redshift surveys are under- 
way. The 2dF Galaxy Redshift Survey (2dFGRS) utilises 
the 2dF instrument at the Anglo- Australian Telescope and 
is based on a revised version of the APM Galaxy Survey 
(Maddox et al. 1990) limited at bj = 19.45. Redshifts have 
now been measured for over 175 000 galaxies (see Colless et 
al. 2001, for a description of this survey). The Sloan Dig- 
ital Sky Survey (SDSS, York et al. 2000) is a CCD imag- 
ing and spectroscopic survey that aims to measure redshifts 
for a sample of 900 000 galaxies. An analysis of the galaxy 
power spectrum from the 2dFGRS is described by Percival 
et al. (2001, hereafter P01). First results on galaxy cluster- 
ing from a subsample of the SDSS are presented by Zehavi 
et al. (2001). 

In addition, a number of other investigations have 
greatly improved the accuracy of various cosmological pa- 
rameters. For example, surveys of high redshift Type la su- 
pernovae have revealed tantalizing evidence for an acceler- 
ating Universe (Perlmutter et al. 1999, Riess et al. 1999); 
the HST Hubble key project has concluded that Ho = 
72 ± 8 km s^Mpc' 1 (Freedman et al. 2001); primor- 
dial nucleosy thesis and deuterium abundance measurements 
from quasar absorption lines imply a baryon density uii, = 
0.020 ± 0.002 (Buries & Tytler 1998ab; Buries, Nollett & 
Turner 2001). With these and many other ambitious projects 
at various stages of development (e.g. weak shear lensing sur- 
veys, CMB interferometers, CMB polarization experiments, 
the MAP, Planck and SNAP satellites^]) it is clear that the 
era of quantitative cosmology has arrived. 

In this paper, we perform a combined likelihood analy- 
sis of the CMB anisotropy data and of the 2dFGRS galaxy 
power spectrum measured by P01. We assume that the ini- 
tial fluctuations were Gaussian, adiabatic and described by 
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power-law fluctuation spectra. Matter is assumed to consist 
of baryons and cold dark matter (CDM) and neutrinos are 
assumed to have negligible rest masses. We allow tensor and 
scalar modes and place no constraints on their respective 
spectral indices and relative amplitudes. Almost all previ- 
ous analyses of the CMB anisotropics have neglected tensor 
modes. However, including tensor modes introduces a ma- 
jor new degeneracy (referred to as the tensor degeneracy in 
this paper) that significantly widens the range of allowed 
parameters (see Efstathiou & Bond 1999, Wang et al. 2001, 
Efstathiou 2001). The tensor degeneracy can be broken by 
invoking additional data sets. Wang et al. 2001 combine the 
CMB data with measurements of the galaxy power spectrum 
from the IRAS PSCz survey (Hamilton, Tegmark & Pad- 
manabahn 2000), estimates of the power spectrum on small 
scales from observations of the Lya forest (Croft et al. 2001) 
and limits on the Hubble constant from HST Hubble Key 
Project. Here we investigate how the major parameter de- 
generacies can be broken by combining the CMB data with 
the 2dFGRS power spectrum. The 2dFGRS power spectrum 
is based on a large survey, with well controlled errors, and as 
demonstrated by P01 already sets interesting limits on the 
matter content of the Universe. Our expectation (see Efs- 
tathiou 2001) is that a joint analysis of the CMB and 2dF- 
GRS will produce accurate estimates of the baryonic and 
matter densities of the Universe and set useful limits on a 
cosmological constant. This expectation is borne out by the 
results described in the rest of this paper. 



2 LIKELIHOOD ANALYSIS 

2.1 Analysis of the 2dFGRS power spectrum 

We use the estimates of the galaxy power spectrum and as- 
sociated covariance matrix computed by P01. As in P01, we 
fit these estimates to theoretical models of the linear matter 
power spectrum of CDM models using the fitting formulae 
of Eisenstein & Hu (1998). The fits are restricted to the 
wavenumber range 0.02 < k/(h Mpc ) < 0.15. Redshift- 
space distortions (see Peacock et al. 2001) and non-linear 
evolution of the power spectrum have negligible effect on 
the shape of the power spectrum at these wavenumbers. 
We will assume that the galaxy power spectrum within this 
wavenumber range is directly proportional to the linear mat- 
ter power spectrum. This is a key assumption in the analysis 
presented in this paper. The lower wavenumber limit is im- 
posed (conservatively) to reduce the sensitivity of the anal- 
ysis to fits to the redshift distribution of galaxies, which are 
computed independently for different zones of the survey. 
Since the 2dFGRS has a complex geometry, the theoretical 
power spectra must be convolved with the spherical average 
over wavenumber of the survey 'window function'. These 
convolved theoretical estimates are used together with the 
spherically averaged estimates of the power spectrum of the 
data and the covariance matrix (computed from Gaussian 
realizations of the 2dFGRS) to form a likelihood function. 
We refer the reader to P01 for a full discussion of each of 
these steps in the analysis. 

In general, the linear power spectrum with wavenum- 
ber measured in inverse Mpc depends on the baryonic and 
CDM physical densities (u;, and lo c ), the scalar spectral in- 
dex n s and an overall amplitude A (the amplitude is treated 
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Figure 1. Contours (1, 2 and 3<r) of the pseudo-marginalizcd 
likelihood functions (see text for details) for various pairs of pa- 
rameters computed by fitting to the galaxy power spectrum of the 
2dFGRS. These contours correspond to changes in the likelihood 
of 2Aln(£) of 2.3, 6.0 and 9.2. The crosses show the position of 
maximum likelihood. 

as an 'ignorable' parameter in this paper and so its precise 
definition is unimportant). However since we use redshift 
to measure distances, the wavenumber of the observations 
scales as h Mpc -1 . The comparison of theory with obser- 
vations therefore requires the introduction of the parame- 
ter h. In fact, the set of variables A, n s , Q m h, u)},/u) m and 
h are natural variables for an analysis of large-scale struc- 
ture: the combination Q m h defines the overall shape of the 
CDM transfer function ( and for negligible baryon density 
is sometimes denoted by the shape parameter V), while the 
ratio uj b /Lj m determines the amplitude of baryonic oscilla- 
tory features in the transfer function (Eisenstein & Hu 1998; 
Meiksin, Peacock & White 1999). 

Fig. 1 shows various two-dimensional projections of the 
'pseudo-marginalized' 2dFGRS likelihood function. When 
using a large number of parameters (as in the CMB and 
CMB+2dFGRS analyses described in the next two subsec- 
tions), it is impractical to compute marginalized likelihood 
contours by numerically integrating over the likelihood dis- 
tribution. Instead, a 'pseudo-marginalized' likelihood func- 
tion in p out of M parameters is computed by setting the 
remaining M — p parameters at the values which maximise 
the likelihood. For a multivariate Gaussian distribution, this 
is equivalent to integrating over the M — p parameters as- 
suming uniform prior distributions (see Tegmark, Zaldar- 
riaga & Hamilton 2001). However, the actual likelihood dis- 
tributions are not exactly Gaussian (as is evident from the 
asymmetrical contours in Figs 1 and 3) and so confidence 
limits assigned to pseudo-marginalized distributions are ap- 
proximate. The contours in the (uj b /u! m , Q m h) plane can be 
compared with Fig. 5 of P01 where the spectral index was as- 
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Figure 2. The points show band-averaged observational esti- 
mates of the CMB power spectrum from Wang et al. (2001) to- 
gether with ±lcr errors. The lines shows the CMB power spectra 
for the adiabatic fiducial inflationary models that provide the best 
fit to the CMB and 2dFGRS power spectra. The parameters of 
these model are listed in Table 1. The solid line shows the best 
fit without a tensor component (fit B). The dashed line shows 
the best fit (fit C) including a tensor component (shown by the 
dotted line). 



sumed to be scale invariant. Relaxing the constraint on the 
spectral index clearly widens the allowed range of LOb/com, 
but the data still place a tight constraint on the 'shape' pa- 
rameter Q m h. As we will see below, the constraints on Q m h 
and n s prove particularly important in breaking degenera- 
cies among parameters inherent in the analysis of CMB data. 



2.2 Analysis of the CMB anisotropies 

The likelihood analysis presented here uses the compilation 
of band power estimates AT§ and their covariance matrix 
Cbb' (including a a model for calibration and beam errors) 
computed by Wang et al. (2001) from 105 CMB anisotropy 
measurements. Each band power estimate is related to the 
power spectrum Ce of the CMB anisotropies by 

ATI = J J2e(i + l)CeW B (e) (1) 

£ 

where Wb is the window function for each band power com- 
puted by Wang et al. These band-power estimates are plot- 
ted in Fig. 2. 

The likelihood analysis of the CMB data uses nine pa- 
rameters. These are: u>b and u) c ; Qa and Qk', the scalar and 
tensor spectral indices n s and n< ; the optical depth to Thom- 
son scattering r opt , assuming that the inter-galactic medium 
was abruptly reionized some time after recombination; the 
amplitude Q of the scalar component and the ratio of r 
of the tensor to scalar amplitudes. Note that definitions of 
the scalar and tensor amplitudes differ from paper to paper. 
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(a) CMB alone (b) CMB alone 
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Figure 3. Results of the nine parameter likelihood analysis. Figs 3a and 3b show approximate 1, 2 and 3ct likelihood contours for various 
parameter pair combinations computed from an analysis of the CMB data alone. Figs 3a use variables natural to the CMB analysis and 
illustrate the geometrical and tensor degeneracies. Figs 3b use the variables natural to the analysis of the galaxy power spectrum (as 
used in Fig. 1). Figs 3c and 3d show the likelihood contours of CMB and 2dFGRS data combined. The crosses in each panel show the 
position of the maximum likelihood. 
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Here we scale the scalar and tensor spectra so that 



x 1(T 5 ) 2 , 



1=2 
50 

±-J2(2t + l)Cl = (2xl0- 5 ) 2 , 



(2a) 
(2b) 



and fit to the data by scaling with the parameters Q and r, 
Ce = Cf + Cj = Q 2 (Cf + rCj). The numbers in equation 
(2) were chosen so that models with Q of approximately 
unity match the data points plotted in Fig. 2 and models 
with f«l have scalar and tensor modes of comparable am- 
plitude. We normalize the spectra in this way to reduce the 
sensitivity of the normalization parameters to other param- 
eters that affect the low order multipole moments {e.g. Oa 
and ilfe) and to decouple Q from the optical depth parameter 
T opt ■ This method of normalizing helps to stabilize searches 
for global maxima of the likelihood functions. For our best 
fit models of Table 1 below we list values of the more com- 
monly used parameter no = Cf /Cf in addition to r. In 
simple models of inflation, the parameters no (or r) , n s and 
tit are related to each other (see e.g. Hoffman & Turner 2001 
for a recent discussion). The relations are model dependent, 
however, and can be violated in multi-field inflation mod- 
els and in superstring inspired models such as the pre-big 
bang (Buonanno, Damour & Veneziano 1999) and ekpyrotic 
scenarios (Khoury, Ovrut, Steinhardt & Turok 2001). We 
therefore assume no relations between no , n s and n t in this 
paper. 

Results from the likehood analysis of the CMB data are 
illustrated in Fig. 3. Almost all of the variance of the param- 
eters used in this analysis, with the exception of Q, comes 
from two major degeneracies (see Efstathiou 2001 for a de- 
tailed discussion) . These two degeneracies are illustrated by 
the likelihood contours plotted in Fig. 3a. The top two pan- 
els illustrate the 'geometrical' degeneracy. This degeneracy 
arises because models with identical matter content, pri- 
mordial power spectra and angular diameter distance to the 
last scattering surface produce almost identical CMB power 
spectra. This leads to a strong degeneracy between Qa and 
Qk, which is broken only for extreme values of Qa by the 
integrated Sachs- Wolfe effect which modifies the shape of 
the CMB power spectrum at low multipoles (see Efstathiou 
& Bond 1999). Since the Hubble constant is fixed by the 
constraint equation, 



h 



(Wfc + 1 



\1/ 2 



(1 - fife - On.) 1 / 2 ' 



(3) 



it is almost unconstrained by the CMB data. 

The lower two panels in Fig. 3a show the constraints on 
the parameter combinations w c -ujb and n s -u>b- These panels 
illustrate the tensor degeneracy: including a tensor compo- 
nent significantly broadens the allowed ranges of parameters. 
For example, values of uib that are more than twice the value 
favoured from primordial nucleosynthesis are allowed by the 
CMB data (Efstathiou 2001). 

Fig. 3b shows likelihood contours using the CMB data 
alone, but computed with the natural variables of the galaxy 
power spectrum analysis as in Fig. 1. The parameter com- 
bination Q m h that essentially fixes the shape of the matter 
power spectrum is extremely unnatural for an analysis of the 



CMB anisotropies. Since Sl m h = (u>b + u) c )/h, the indeter- 
minacy in h arising from the geometrical degeneracy smears 
the likelihoods along the direction of fi m /i. The wide range 
of allowed values of tub /u m and the tight correlation with n s 
is a consequence of the tensor degeneracy. 

2.3 Combining the CMB and 2dFGRS likelihoods 

Fig. 3b is interesting because it shows that the CMB like- 
lihoods in three of these plots are complementary to those 
of the 2dFGRS analysis (LOb/ui m -£l m h, n a -Q, m h, LUb/^m - 
n s ). The addition of the 2dFGRS constraints breaks both 
the geometrical and tensor degeneracies, resulting in strong 
constraints on u>b, uj c , ^a and h. The way that this works is 
evident from Figs 1 and 3b: the constraints on n s from the 
2dFGRS help to break the tensor degeneracy by excluding 
high values of uib and low values of lo c . The resulting values 
of uJb and u> c fix the Hubble radius at the time that mat- 
ter and radiation have equal density, which in turn largely 
fixes the shape of the CDM transfer function in physical 
Mpc. Comparing with the power spectrum of the 2dFGRS 
in 7i _1 Mpc constrains the Hubble constant, thus breaking 
the geometrical degeneracy. 

The lower panels in Fig. 3 show the results of combining 
the CMB and 2dFGRS likelihoods. The results are striking, 
showing a significant tightening of the constraints in each 
plot. Table 1 lists parameters corresponding to maximum 
likelihood fits to the data and the approximate ±2<r ranges 
of each parameter. The second column lists the maximum 
likelihood fit to the CMB alone (fit A). The parameters of 
this fit are identical whether or not a tensor component is 
included. The third and fourth columns (fits B and C) list 
the maximum likelihood fits to the CMB and 2dFGRS data 
excluding and including a tensor mode. The fifth column (fit 
D) adds the constraint from big-bang nucleosynthesis (BBN) 
of a Gaussian distribution for u>b centred at u>b = 0.020 with 
a dispersion of Au>b = 0.001 (Buries et al. 2001). 

The parameters of fit B, which provides a perfectly ac- 
ceptable fit to the data, are very close to those of the stan- 
dard 'concordance' cosmology (e.g. Bahcall et al. 1999). In 
particular, the baryon density is compatible with the pri- 
mordial nucleosynthesis value, and the Hubble and cosmo- 
logical constants are compatible with more direct observa- 
tional estimates. The CMB power spectrum for this solution 
is plotted as the solid line in Fig. 2 and the linear matter 
power spectrum is plotted together with the 2dFGRS data 
points in Fig. 4. Both curves provide acceptable fits to the 
data. Fit B has a low baryon fraction of LO b /uJm = 0.15. As a 
consequence, the amplitudes of the baryonic features in the 
matter power spectrum are almost imperceptibly small (see 
Fig. 4). 

Allowing a tensor component produces a slightly better 
fit to the data, but the parameters are less concordant with 
other observations (Fit C, Table 1). The CMB power spec- 
trum for this model is plotted as the dashed line in Fig. 2. 
According to this solution, a significant part of the COBE 
anisotropies comes from a tensor component. The baryon 
density of fit C is co b = 0.027 and is well outside the range 
of values inferred from primordial nucleosynthesis. The mat- 
ter power spectrum for this model is plotted as the dashed 
line in Fig. 4. This shows clearly what is happening with 
this solution. The apparent wiggles in the 2dFGRS power 
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Table 1: Parameters values and errors 



Fit A 
CMB alone 
+ tensor 



Fit B 
CMB+2dFGRS 
no tensor 



Fit C 
CMB+2dFGRS 
+ tensor 



Fit D 

CMB+2dFGRS+BBN 
+ tensor 



approximate ±2cr paramctc 

Fit A Fit C 

CMB alone CMB+2dFGRS 
+ tensor + tensor 



ranges 
Fit D 

CMB+2dFGRS+BBN 
+ tensor 



T~c.pt 

nt 

r 
h 



0.020 
0.13 
0.96 

-0.04 
0.43 




0.14 



0.021 
0.12 
1.00 
0.001 
0.71 




0.15 
0.21 
0.69 



0.027 
0.085 
1.20 
-0.030 
0.80 


-0.10 
0.60 
1.24 
0.24 
0.16 
0.71 



0.020 
0.10 
1.04 
-0.013 
0.73 


0.13 
0.20 
0.26 
0.17 
0.19 
0.66 



0.016-0.045 
0.03-0.18 
0.89-1.49 

-0.68-0.06 

< 0.88 

< 0.5 

< 0.98 



0.018-0.034 
0.07-0.13 
0.95-1.31 

-0.05-0.04 
0.65-0.85 

< 0.5 

< 0.87 

0.13-0.28 
0.12-0.22 
0.60-0.86 



0.018-0.022 
0.08-0.13 
0.95-1.16 

-0.05-0.04 
0.65-0.80 

< 0.5 

< 0.82 

0.13-0.22 
0.16-0.21 
0.61-0.84 
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Figure 4. The points show the galaxy power spectrum of the 
2dFGRS measured by P01 divided by the power spectrum of a 
scale-invariant CDM model with £jj = 0, Q m h = 0.2. The er- 
ror bars are computed from the diagonal components of the co- 
variance matrix. The lines show the linear matter power spectra 
of the maximum likelihood fits to the combined CMB and 2dF- 
GRS power spectra after convolution with the spherically aver- 
aged window function of the survey. The solid line shows fit B 
from Table 1 (no tensor component). The dashed line shows fit C 
(including a tensor component). 



spectrum pull the solution towards a high baryon fraction. 
However, to produce a good fit to the CMB anisotropies 
with a high baryon fraction, the tensor degeneracy of Fig. 3 
requires high values of n s and significant tensor anisotropies. 
The likelihood ratio of fits B and C is Lb / £c = 0.34 and 
so fit C is only marginally preferred over fit B. In two of the 
panels from Fig. 3c and d, the likelihood distributions have 
two peaks centred at the parameters of fits B and C. Adding 
the BBN constraint on Ub (fit D) selects one of these peaks 
with parameters close to those of fit B. 

Fits B and C predict a lower normalization for the 
present day matter power spectrum than implied by the lo- 
cal abundance of rich clusters of galaxies. In a recent anal- 



ysis of the number density distribution of rich clusters as a 
function of X-ray temperature, Pierpaoli, Scott and White 
(2001) deduce 



<T8 



(0.495 



+0.034 v n - 
0.037 j s 'n 



(4) 



where erg is the rms fluctuation in the mass density distri- 
bution averaged in spheres of radius 8ft Mpc. Fit B gives 
erg = 0.72 and fit C gives erg = 0.61, whereas equation (Q) 
implies erg = 1.04 and erg = 1.20 respectively. Most of the 
error in equation (Q) comes from uncertainties in the cluster 
mass- X-ray temperature relation and it is not clear whether 
the quoted error reflects the true uncertainties. A number of 
effects could boost the best fitting values of erg, for example, 
a realistic value for r opt or possible calibration errors in the 
CMB data might affect erg at the 10-20% level. Such effects 
may reconcile fit B with the cluster data, but are probably 
not large enough to explain the discrepancy with fit C. Fur- 
thermore, as we have discussed above, the discrepancy with 
the primordial nucleosynthesis value of provides another 
reason to disfavour fit C. 



3 DISCUSSION 

The results of this paper are based on the key assumption 
that the galaxy power spectrum on large scales (wavenum- 
bers k < 0.15ft Mpc" 1 ) is proportional to the linear matter 
power spectrum. Under this assumption, we have shown that 
the galaxy power spectrum of the 2dFGRS can be used to 
partially break the two major parameter degeneracies inher- 
ent in the analysis of CMB anisotropies. The limits on the 
scalar spectral index from the 2dFGRS help to break the ten- 
sor degeneracy. The resulting constraints on the matter den- 
sity provide a measure of a standard physical distance (the 
Hubble radius at the time that matter and radiation have 
equal density). This standard length constrains the Hubble 
constant and so breaks the geometrical degeneracy. 

The resulting constraints are in remarkable agreement 
with the baryon density inferred from primordial nucleosyn- 
thesis (Buries and Tytler 1998ab), estimates of the Hubble 
constant from the HST Hubble key project (Freedman et 
al. 2001) and evidence for a non-zero cosmological constant 
from obserations of distant Type la supernovae (Perlmutter 
et al. 1999; Riess et al. 1999). The best fit model excluding 
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a tensor component has parameters that are very close to 
those of the standard 'concordance' cosmology (Bahcall et 
al. 1999). However, the combined CMB+2dFGRS data pro- 
vide weak upper limits on a tensor component (Table 1) and 



other so lutions arc possible which have a higher baryon den — 
sity an a baryon fraction. Thoao colutionc conflict with tho 

limits on uib from primordial nucleosynthesis and require a 
scalar spectral index n s > 1. The model with high n s and 
high u)b provides a somewhat closer match to the appar- 
ent 'wiggles' in the galaxy power spectrum at wavenumbers 
k ~ O.O8/1 Mpc -1 and k ~ 0.12ft Mpc -1 than is achieved by 
the scalar only model (Fig. 4). Neither model fully matches 
the data points, however, and it is plausible that the ap- 
parent features are enhanced by the noise. New power spec- 
trum data from the 2dFGRS and the SDSS will soon al- 
low us to test this hypothesis. It is particularly encourag- 
ing that the combination of the 2dFGRS and CMB data 
yields strong evidence for a cosmological constant in the 
range 0.65 ;$ Qa ,$,0.85 based on completely different argu- 
ments to those applied to distant Type la supernovae. This 
significantly strengthens the case in favour of an accelerating 
universe. 
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